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THE ROLE OF DROP VELOCITY IN STATISTICAL SPRAY DESCRIPTION 


J. F. Groanawag,^ M. N. Bl*Hakll,^ F. 8. Myari,^ and 0. A. Uyahara^ 


ABS TRACT 

The Juscificaclon for describing a spray by 
creating drop velocity as a random variable on an equal 
statistical basis with drop size was studied experimen- 
tally. A double-exposure technique using fluorescent 
drop photography was used Co make size and velocity 
measurements at selected locations in a steady ethanol 
spray formed by a swirl atomizer. The size-velocity 
data were categorized Co construct bivariate spray 
density functions to describe Che spray inxnediaCely 
after formation and during downstream propagation. It 
was found chat a statistical treatment of drop velocity 
was supported by Che data. Spray density function 
shapes and modal characteristics depended strongly on 
position and the amount of droplet-gas interaction that 
had occurred. Bimodal density functions were formed 
by environmental interaction during downstream propaga- 
tion. Large differences were also found between spa- 
tial mass density and mass flux size distributions at 
the same location. 

NOMENCUTTRE 

D ■ drop diameter, am 

f • spray density function, no. drops per unit 
size, velocity, position, drop temperature 
fg * measured spray density function, no. drops/ 

(um- cm/s ' cm^ ) , Eqs . (11) and (12) 

■ flux drop size distribution, no. of drops/ 

(am ■ cm^ • s ) 

fg • spatial drop-size distribution, no. of drops/ 
(um- cm^ ) 

fy * velocity distribution, no. of drops/(cm/s) 

■ onc-dlmenslonal spray density, no. of drops/ 
(am-cm-cm/s) , Eqs. (12) and (18) 

g • function of drop variables, Eq . (2) 

M • drop mass, Cj^irD^/6, g/ cc 

m ■ constant In Eq. (13) 

N ■ total no. of drops 

"ijk “ drops in 1^^ size, axial velocity, 

and radial velocity categories 

r ■ radial coordinate, cm 

Tj_ ■ temperature of drop, K 

C ■ time, s 

^ Lewis Research Center, National Aeronautics and 
Space Administration, Cleveland, Ohio 
^ University of Wisconsin, Madison, Wisconsin 


u - gas velocity in spray, cm/s 

V, ■ sampling volume, cuP 

Vg - liquid exit velocity, cm/a 

Vjj • axial velocity of liquid sheet, cm/s 
V ■ drop velocity vector, cm/s, components in 
cylindrical coordinates, Vj., v^, v^ 
w > weighting function, Eq. (2) 

X - position vector, cm 

z ■ axial coordinate, ran 

e > no. of photographic samples at given spray con- 

dition 

• general drop variable 
at • Increment in 

■ drop density, g/cm^ 

“ spray density, g/cm^ 

_ - underline -- vector 

^ / « ensemble average 

( • ensemble average weighted by M 

- average in category 
' • particular location 

INTRODUCTION 

In the hundreds of spray studies which have been 
conducted over the past century droplet size has re- 
ceived the most emphasis as the key variable in spray 
description. It has long been recognized that the 
atomization process Is random as far as the sizes of 
the droplets formed and a statistical treatment in 
terms of size distributions and associated means is 
used. The situation with respect to other droplet 
Independent variables such as velocity is much less 
clear. 

A general statistical mechanical theory of sprays 
has been formulated [ij which includes size, velocity, 
and position as Independent variables, but applications 
of the theory have been primarily limited to consider- 
ations of the statistical properties of drop size. 

While some measurements of drop velocity have been 
made [2, 3,4,5], velocity has not been purp sely 
treated as a random variable on an equal statistical 
basis with drop size. Data have been interpreted on 
the supposition that drops of a given size all move 
with the same average velocity at a given location In 
the spray. 

However, physical examination of the spray situa- 
tion leads to the conclusions Chat (a) spray formation 
la a random process which Is distributed In space and 
(b) the histories of Individual droplets are unique 
functions of the initial conditions and later environ- 
ment. These conditions Imply that, at downstream lo- 
cations wtiere a collection of approximately spherical 
drops first exists, other droplet properties besides 
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■Ize should be created as being scaclaclcally dlstrib- 
uCed. 

'"'e Inves clgatlon described In this paper ex- 
ploi > .he Justification for and consequences of treat- 
ing drop velocity as a random variable on an equal sta- 
tistical basis with drop size. The approach was basic- 
ally an experimental one In which detailed measurements 
of drop size and velocity were made at selected loca- 
tions In a spray. The measurements were snalyzed In 
terms of a spray density function which Is based on a 
statistical mechanics approach to spray description. 
Specific goals of the study were to determine xf dsca 
support such a statistical treatment of drop velocity 
and to examine the Implications of such a treatment 
for Interpreting size data and analyzing spray behav- 
ior. Because of the difficulties In measuring sizes 
and velocities and the large number of measurements 
required, Che scope of the Investigation was limited 
CO a detailed characterization of one particular spray 
situation. 






hNf dD dv dT^ 


Mf do dv dT^ 


where ^vj^ remains a function of position x ~ x* 
and M la the droplet mass, PjTtiP/ 6. The denominator 
■ Ps = 


of Eq. (3) la Che spray density p, 

/// Mf dD dv dTL 




(4) 


Various marginal density functions may be defined 
In which all Independent variables except one are In- 
tegrated out mathematically or disregarded experimen- 
tally. The marginal density on la: 


f(Pi) -yf(r^,rj,t)di'j 1 .1 J (5) 

where, as usual, the Integrals are over the whole 
range of Pi. A particular marginal density Is Che 
usual spatial drop size distribution: 


SUlWARY OF STATISTICAL SPRAY DESCRIPTION 

For purposes of this discussion, spray studies 
are divided Into the two general areas of formation 
and propagation. Formation concerns Che process by 
vdilch liquid In a reservoir Is atomized. The descrip- 
tion of Che resultant spray depends upon Initial con- 
ditions of fluid properties, atomizer geometry, energy 
addition and properties of the medium In which Che 
spray Is formed. Propagation Involves Che description 
of changes In spray properties due to transfer pro- 
cesses which occur In the two-phase flow downstream. 
The conceptual boundary separating these two regimes 
may be called the surface of formation. Describing 
the spray at this interface Is the endpoint of forma- 
tion studies and the Initial condition for propagatloi. 
studies. The method of describing Che spray In terms 
of a spray density function at the surface of forma- 
tion and some aspects of Its propagation are the main 
emphases In this paper. 

Spray Density Function 

Consider a function f(r.,c) as a representation 
of the spray at Che surface of formation. The varia- 
bles are randomly distributed and are chosen as 

chose properties necessary to describe the droplet 
state at formation and during subsequent propagation, 
that Is, drop mass, momentum and energy at any loca- 
tion. Therefore, reasonable choices for are size 

D, position X, velocity v, and temperature Tl- The 
Inclusion of Tj, Is a generalization of the treatment 
In reference 1 and Is discussed In detail In refer- 
ence b. The function f represents the probable num- 
ber of drops In the range dP^ about at a time 

t. In this unnormalized form: 



where N Is the total number of drops represented by 
f at a time t. 

Ensemble averages or moments of a function 
g(r^,t) when weighted by another function U(r^) are 
given by: 

.t);„ - 


For example, the :.dss average velocity is given by: 


f w(rt)g(rj.t)f(rj,t)dri 

f wcr^jfcrj.tjdi'i 


( 2 ) 


fg(D,x',t) “ JJ t(D.v,x' .T^.tjdv dTj^ (6) 

Similarly, the "velocity distribution" Is: 

fv(v.x’,t) ■// f(D.v,x'.TL,t)dD dTL 

Another marginal size distribution of Interest Is the 
temporal or flux size distribution given by: 

fF(D>i'.t) - JJ vf (D.x* ,v,TL,t)dv dTL 

The physical significance of fg and fp are as fol- 
lows. The spatial distribution Is the number of drops 
per unit spstlal volume per unit size at a time t 
and Is measured by Instantaneous photographs of drops 
In a known spatial volume. The flux distribution Is 
the number of drops per unit area per unit size and 
time collected as they cross a known surface area 
averaged over time. Equation (8) defines the flux 
distribution In a way which Is equivalent to the one 
obtained by collection If the spray Is a stationary 
random process with time averages equal to ensemble 
averages . 

The average (expected value) of any function 
g(r.) for a given value of the spray variable I • Is 
defined for 1 J as; 

- 

A particular case of Interest Is the expected value of 
velocity at a given .ilze as a function of D and x' 
given by : 

lilo; ■ 


^^/~vf(D|V|X',TL.t)dv dTL ^F 

yy’f(D,v,x',TL,t)dv dTL 


(10) 


7g(rL)f(rprj.t)drL 

J f(ri,rj,t)dPL 


(9) 


It has been estimated experimentally by averaging mea- 
sured velocities of particular sized drops [3,4,3] but 
was not Interpreted ss being derived from s density 
function f which contained v as a randomly distrib- 
uted variable. Equation (10) also shows that (v|0) 

Is given by the ratio of flux to spatial size distri- 
butions . 

This brief stimmarv of density functions, derived 
forms and means Is sul lie lent for the Interpretation 
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of the data to be presented. For a more detailed die* 
cuasion of equations of change for f and aaaoclacad 
owana, aee Refs, [l] and [6j. 

MEASUREMENT OF THE SPRAY DENSITY FUNCTION 

The particular form of the denalty function mea- 
sured was restricted to the variables D, y, and x. 
Droplet temperature was not measured which was equiva- 
lent to Integrating over T^. Injection parameters 
were held constant so that a steady-state spray condi- 
tion was assumed eliminating consideration of time, t. 
In terms of the general f, the measured density fg 
Is: 

fg(D,y,x') - y f (D,y,x' ,TL,t)dTL (U) 

where x' Indicates that the functions are evaluated 
at a particular spray position, x > x'. Synmetry was 
assumed about the spray axis, and only two components 
of position and velocity were considered In cylindri- 
cal coordinates: r, z, v^., and v^ - Any v^ was 

normal to the sampling plane and could not be mea- 
sured . 

Fluorescent Spray Photography 

A double-exposure fluorescent technique was used 
to measure sizes and velocities. This method left the 
spray undisturbed and provided direct local values of 
the droplet variables. The fluorescent technique of 
photographing droplets was originally developed [7,8] 
as a single exposure method of measuring the sizes of 
drops In a small spatial volume at any Instant to give 
the spatial drop size distribution, f . In the pres- 
ent study a precisely controlled double-exposure capa- 
bility was added to provide a measure of droplet ve- 
locity as well as size so that fg could be esti- 
mated . 

The key feature of the technique Is the addition 
of a fluorescent dye to the liquid being sprayed. A 
shaped light beam selectively lights a region of the 
spray and defines the sampling volume on which the 
camera Is focused. Only drops within the camera's 
depth of field are lighted and caused to fluoresce. 
Multiple exposures may be recorded without loss of 
contrast since each fluorescing droplet Is a primary 
source recording Its Image with an unllghted back- 
ground . 

A pictorial view of the experimental arrangement 
which was used Is shown In Fig. 1. The axes of the 
spray, camera and lighting system are mutually perpen- 
dicular. Liquid containing the fluorescent dye was 
Injected vertically from a swirl atomizer, passed 
through the region where sampling occurred, and was 
collected and removed from the room by an exhaust sys- 
tem Identical lighting systems consisting of con- 
stricted spark gaps and quartz condensing lenses hav- 
ing specially shaped aperture stops were located on 
either side of the camera axis. When either gap was 
fired the condenser lenses focused and shaped the beam 
to light the same volume In the spray. The firing 
sequence of the two sources was monitored by a photo- 
tube and controlled to produce two flashes separated 
by a known time Interval. Each drop within the sam- 
pling volume viewed by the camera succesalvely fluo- 
resced and was recorded twice on Che film. The posi- 
tion of the camera and lighting system was fixed to 
maintain alignment, and the nozzle was positioned so 
that the spray could be sampled at various axial and 
radial locations. 

Table I summarizes the specifications and oper- 


ating conditions for the various elements of the sam- 
pling system. Details of Che development of the light- 
ing and camera systems and Che fluorescent dye charac- 
teristics are available elaawhara [7,8]. Successful 
application of the fluorescent technique depends on 
having very Intense light sources and carefully align- 
ing the system elements. 

Sampling Conditions and Pats Reduction 

An unconflned sprav was formed by steady Injection 
of ethyl alcohol through a swirl scomlzer Into s room 
at ambient temperature and pressure. The sampling lo- 
cation and Injection pressure were varied as shown In 
Che diagram of sampling geometry given In Fig. 2. A 
traverse of radial positions at locations linnedlacely 
after breakup showed Initial conditions Just after 
spray formation for three Injection pressures. Radial 
surveys at two downstream distances provided data on 
the changes thec occurred during propagation. The 
magnitude of Che air velocity In Che spray was due to 
a combination of motion Induced by the exhaust fan, 
which was less than 60 cm/aec over Che range of posi- 
tions measured, and the entrained air motion produced 
by the momentum transfer from Che spray to the air. 

The photographs were taken with the room darkened 
and several samples of collections of droplet Image 
pairs were recorded on each film. The Interval be- 
tween flashes was chosen so chat complete separation 
of Che Images was achieved for most of the Image pairs. 
A diameter, a separation distance and an angle of Che 
drop trajectory in Che r-z plane were measured for 
each droplet pair. The measurements were made by hand 
on the screen of a microcard reader at addlclo.ial mag- 
nifications of 16 to 18 times. From these measure- 
ments the size D and two components of velocity, Vf 
and Vj, were calculated for each drop. The conical 
liquid sheet velocity before breakup was also meas- 
ured from double exposure photographs. A total of 
nx>re chan 30,000 pairs of drop Images were measured 
and the data were processed on a digital computer. 


ANALYSIS OF THE SIZE-VELOCITY DATA 


The dem Ity function fg is estimated by cate- 
gorizing the size-velocity data and applying an ap- 
proximate form CO give values of f^; at Che category 
moans : 


*E^®l'''zJ '’rk 


,x') 






'rk 


(12) 


The estimate Is Improved as Che sample size 
the number of samples e Increase and the category 
sizes are decreased. Category boundaries were chosen 
which Increased by a constant multiple so that the 
fractional change in D or v was a constant. For 
any drop variable: 


Pi - V-l 


(13) 


T, 


l+l 


( 14 ) 


and the category geometric mean was used: 

Ti • (rjrj^i)^'^ (15) 

The values of m and chosen for D, v^, and v^. 

were, respectively, 1.31, 1.31 and 1.73 for m and 

10 urn, 38.1 cm/s and 123.4 cm/s for the first boundary, 
Pj^. Once values for fj- were obtained, any of the 
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marginal denslcies auch as weighted denalcles 

auch aa or mean quantltlea auch aa can be 

calculated by summing over the categories to approxl* 
mate the Integrals defined previously. 

A one-dlmenslonal density function, f^, was ob> 
talned by Integrating fg over the cross section. In 
cylindrical coordinates with equal radial Increments, 
the Integral is approximated as a sum over I radial 
stations ; 


. V CZ: - 0 n 

Vs -=^1 ^''sj ‘^''rk y 


(lb) 


Radial points were taken at 2 tim Increments beginning 
at r • I mm. 

Categorization of the data gave estlmatea of the 
density functions. In this paper the densities are 
plotted as continuous curves through data points at 
the category means rather than histograms. Normalized 
forms are used weighted by droplet mass, M. The nor- 
nallzed form allows comparisons of density shapes to 
be made at widely different conditions and Che mass 
weighting emphasizes the larger sizes avoiding ex- 
treme skewing Coward small sizes characteristic of 
number densities. 


Marginal Density Functions 


Spatial and Flux Distributions 


Blmodal spatial drop size distributions obtained 
by photographic methoda have been reported by several 
Investigators. The most similar study Co the present 
one used a swirl atomizer Injecting Into stagnant air 
In a closed chamber [llj. Measured values of fg 
were strongly dependent on location, and In many cases 
were decidedly blmodal. Atomization by Impinging Jets 
Injecting Into still ale [12] and higher velocity air- 
screams [14,15] have also produced spatial distribu- 
tions with two modes. Due to Che difficulties In 
separating true modes from statistical fluctuations In 
small samples and the complexity of treating blmodal 
data analytically, much data has been assumed to be 
unlmodsl [8]. It Is probable that reanalysla of much 
existing photographic data would reveal the existence 
of two distinct modes. 

There Is also a body of data obtained by collec- 
tion methods [e.g. 13] or velocity weighting of spatial 
distributions which corresponds to the flux distribu- 
tion, Aa stated In Eq. (10) the ratio of ^ to 

fg Is the average drop velocity at a given size: 

(v|D). Figure 4 compares Che two normalized distribu- 
tions at a particular downstream location. From Eqs. 
(3)> (^)t and (10), the ratio of the mass-weighted, 
normalized forms Is: 


Mfpz /fj_\ _ (Vz’D) 

f's(''z;M\Mfsj‘ (''z-M 


(17) 


The three marginal mass weighted densities for 
D, V 2 , and v^ are shown in Fig. 3 at the lowest In- 
jection pressure near Che surface of formation. Varia- 
tions with radial position are shown for each case. 

The curves have been drawn Co emphasize their essential 
blmodal character with fairing through the scatter at 
large sizes. The spatial density In Fig. 3(a) Is the 
usual spatial drop size distribution obtained by photo- 
graphic sampling. Note Che strong radial dependence 
with the mode at lower sizes Increasing In prominence 
near Che outside of the spray. 

Figure 3(b) shows the same spray condition with 
axial drop velocity as Che Independent variable. As 
before, the development of an unambiguous second mode 
at lower velocity is shown at the outer radial loca- 
tions. As will become clear, this mode corresponds to 
drops being decelerated to local air velocity. Final- 
ly, Fig. 3^c) shows the marginal mass density as a 
function of radial velocity. Two distinct modes are 
again present at each location. The mode at smaller 
velocities peaks at small negative values of v^. re- 
lated to the inward flow of entrained air. Similar 
overall formation behavior was observed at other Injec- 
tion pressures . 

The entire local behavior of f^ Is very diverse 
and spatially dependent, local values of Che bivariate 
forms of fg(D,V 2 ,r ' ,z ' ) could be plotted with para- 
metric cuts through Che surface at constant D or 
constant v^ . Their overall character Is similar to 
Fig. 3. Tlie density properties are Initially con- 
trolled primarily by Che atomizer used and the spray- 
ing parameters. However, the spatial denalcles arc 
radically modified by drop-gas Interactions (mainly 
drag in this case) as propagation proceeds. Local 
values of air velocity determine the location of the 
developing low size-velocity mode which becomes more 
pronounced with Increased travel time from the forma- 
tion region. Ttie coupling of the liquid flow with the 
gas produces .ilr entrainment [9,10]. Thus, relative 
velocity, which Is the driving force for changes in 
fg. Is a function of position. 


as Is shown In Fig. 4. In this case small drops have 
decelerated and their spatial density has Increased 
while Che largest drops continue Co move much faster. 
Thus, for this gas flow condition photographs show 
the largest population of small drops while collectors 
Intercept a greater number of large drops . Figure 4 
emphasizes the fact that the two size distributions 
arc not equivalent and may differ substantially. Only 
In the special case where all drops are traveling at 
the same velocity are the normalized forms of and 

f^ equal. Figure 4 also shows that the modal charac- 
teristics of and fg can be drastically differ- 

ent. At this downstream position the spatial size 
distribution Indicates Chat drops less than about 
60 4m have nearly reached Che air velocity while 
larger drops continue to n.ove faster and account for 
the mode at large sizes. But wlien Che spatial density 
Is weighted by the velocity regression curve to give 
mass flux, Che dominant mode appears at large sizes 
with only a small Inflection remaining In the small 
size range. The collection of regression curves 
(divided by different radii at a down- 

stream location appear In Fig. 5. It can be seen that 
these weighting curves which relate Che two types of 
distributions reflect the stage of deceleration so 
Chat their range varies greatly from the Inner Co 
outer locations In Che spray. 

One- Dimensional Density Functions 

A one-dlmenslonal description of the spray at 
any axial location Is obtained by Integrating f over 
a cross section. In cylindrical coordinates: 



The corresponding numerical approximation using ex- 
perimental data la given by F.q . (lb). Note that f^ 
Is a one--'lmenslonal "density" with units of drops 
per unit size, axial velocity and length In the z 
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direction; and fluxes obtained from f^ are simply 
flow rates In the axial direction. 

Normalized mass densities as a function of D 
are shown in Fig. 6(a) for the location near the sur- 
face of fomuition at the three injection pressures. 

The small first modes show the influence of the outer 
regions of the spray where drop deceleration is ap- 
preciable. A shift toward smaller sizes with increas- 
ing Ap, illustrates the well known fact that the 
higher energy inputs produced smaller drops. Figure 
6(b) shows the corresponding velocity dependence. 

Only small fractions of the mass have approached 
equilibrium with the air. The exact location and 
height of the low velocity mode in each case depends 
on how well the sampling location approximated the 
surface of formation and the extent of the formation 
region which was greatest at .\p - 1.7 atm. Nearly 
all of the mass is located in the second modes which 
broaden with increasing Ap and have means ranging 
from 2/3 to 3/4 of the axial sheet velocity. 

The propagation of the mass densities with down- 
stream distance is traced in Fig. 7. With respect to 
either size (Fig. 7(a)) or velocity (Fig. 7(b)) there 
is the progression from a dominant second mode, 
through modes of comparab,.e size, to a dominant first 
mode as more and more of the mass approaches gas ve- 
locity. 

A clearer picture of the clianges occurring during 
propagation is given in Fig. B) .by the contour plots 
of the bivariate size-velocity function. Inmedlately 
after formation (Fig. W(a)) the large "lilll" repre- 
senting the second, or wliat may be called the forma- 
tion mode, is dominant. Only a small peak representing 
Che first, or propagation mode, appears. At the medium 
downstream distance (Fig. 8(b)) Che propagation peak 
has sharpened; the formation mode has diminished; and 
a higher ridge connects the two. Finally the dominant 
feature of Fig. 8 (c) at Che farthest downstream dis- 
tance Is rhe high propagation mode witose base blends 
into the extensively altered fermatlon mode. 

It is if interest to note Chat in the studies 
wliere blmodal distributions were reported from samples 
taken at a constant downstream distance [14,15] it was 
Che large size node which showed the usual changes at- 
tributable to variations In injection parameters. 

This is compatible with Che concept that it repre- 
sented the formation mode while the first mode indi- 
cated Che stage of propagation. Tlie alternate hypoth- 
esis that the two modes resulted from two distinct 
formation processes is possible, and only velocity 
data could decide the question. 

In the present study the possible existence of 
blmodal formation processes can be investigated by 
considering the one-dimensional mass flux distribu- 
tions. Since vaporization was small, £p, the axial 
flux drop size distribution, should propagate nearly 
unchanged with z. If definite modes are present at 
formation they should appear in the flux distribution. 
Figure 9 , wtilch Includes both formation and propaga- 
tion Information, indicates that no definite modes are 
present for Ap ■ 1.7 atm. Ttie strong blmodal spatial 
characteristic has almost completely disappeared at 
downstream locations. It Is seen that the flux dis- 
tribution docs remain approxlmrtely constant for the 
three locations with the small shift toward smaller 
sizes probably being caused by vaporization. Little 
difference exists between the two curves for higher 
pressures. The position of the two points at 75 and 
100 urn for 2.72 atm could be interpreted as a blmodal 
formation tendency. However, the data are not exten- 
sive enough to warrant a definite conclusion. Wliat Is 
clear Is that the two modes considered throughout the 


discussion of spatial density functions are the result 
of the drop-gas interactions and are not inherent in 
the spray formation process. 

The formation and propagation results from Figs. 
6, 7, and 8 are suaMrlzed schematically in Fig. 10 
for the situation existing in the present study wliere 
the mean velocity at the exit of the atomizer, Vg, ia 
greater than the ambient air velocity, u. Each con- 
tour plot of the mass density surface in the size- 
velocity plane is accompanied by the marginal density 
-functions of 0 or v alone. The relationship be- 
tween the spatial and flux drop size distributions Is 
also noted for each case. At the formation and equi- 
librium conditions the flux and denalty diatrlbutlons 
differ by approximately a constant. Density contours 
show the progression from a single "hill" formation 
mode to the blmodal intermediate propagation surface 
and, finally, the return to a single hill where drops 
and gas are in velocity equilibrium. The ideal equi- 
librium condition where the velocity density approach- 
es a delta function was not yet reached at the actual 
downstream locations sampled (Fig. 8(c)). Reasoning 
from Fig. 10, the corresponding curves for the case 
of injection into a higher velocity gas stream can be 
visualized by Inverting the contours and densities 
with respect to the velocity ordlnacc. Drops would 
then accelerate from Vg to u. 

CONCU'S IONS 

Detailed drop size and velocity measurements 
made on a spray immediately after formation and I'uring 
propagation lead to the following conclusions about 
the rule of drop velocity in spray description: 

1. Drop velocity in a spray Is a statistically 
distributed variable, the knowledge of which is as Im- 
portant as drop size. Beginning at the shortest down- 
stream distances wliere a spray may be said to exist, 
the velocities of drops of a given size are distrib- 
uted over a range of values usually with a mean sub- 
stantially different from the mass average injection 
velocity . 

2. Drop concentrations In particular size and ve- 
locity ranges as given by the size-velocity density 
function are strongly dependent on position. Tills 
spatial dependence is determined to a large extent by 
the amount of droplet-gas Interaction that has oc- 
curred. A wide variation In the shape and modal char- 
acteristics of the distribution curves may be found In 
the same soray simply by sampling at ditferent loca- 
tions. I'nless sampling conditions are very carefully 
specified, little basis exists lor the comparison of 
ditferent sets ot spatial drop size data. 

3. Blnuidal density functions may be formed by en- 
vironmental interactions occurring during downstream 
propagation. In the present case, tie selective de- 
celeration of drops according to size produced a tran- 
sition from a primarily unlmodal formation character- 
istic to a blmodal condition consisting of drops near 
velocity equilibrium with the gas and Chose still be- 
ing decelerated. 

4. In many cases, large ditlcrences in shape and 
modal characteristics exist between spatial and flux 
distributions at the same location. Unly In special 
situations where all drops have closely approached the 
same velocity are photographic and collection data 
equivalent. 
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TABU 1. - SAMPLING SYSTEM SPECIFICATIONS AND 
OPERATING CONDITIONS 

Llaht aourcaa ; Guided air aparka In 0.8 x 19 nn slot 
Maximum anargy • 80 Joulat; 0.1 uf charged to 40 kV 
Flaah duration (half peak) • 1.5 to 2.0 ua 
Dalay between flaahes - continuously variable; 

nominal valuaa uead 9.5 to 74 us 
Control - spark gap twitches pressurized with dry 
nitrogen 

Sampling volume ; Thin slab parallel to the spray axis 
Size - 4.0 X 4.0 x 0.2 mn 

Formed by - two 152 nzn f/1.2 plano-parabollc, 
fuaed-quartz condensing lensea 
Fluorescent dye ; Uranln (fluorescein) 

Concentration - 5 grama/llter In 95% ethyl alcohol 
Spactral characterlstlca - absorption peaks at 
2500 and 4900 X; emission peak at 5300 X. 

Camera ; Two lent ayttem 

Objective lens - f/3.5, 152 nin operated at f/5.6 
Relmaglng lent - f/2.0, 58 nin 
Overall magnification • 25 

Size resolution - lOu ± 10% (static calibration) 
Working distance - 152 inn 
Depth of field for lOu objects - *>220u 
Film : 4x5 sheets, ASA 1200 

Development - 12 minutes continuous agitation 
followed by Monckhoven's Intenalfler 
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Figure 1. - Experimental arrangement for double-exposure fluorescent photography. 
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Figure 2. - Sampling geometry and conditions. 
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Figure \ - Mass densities near the surface of formation. 
1.7 atm. 7 • 15.9 mm. 
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Figure 4. - Comparison of the modal characteristics of 
spatial and flux dis^rtouticns at a particular location. 
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gure 5. - Expected values of drop velocity as a function of size at sev- 
eral radii at a downstream location. 
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Figure 6. ■ One-dimensional mass densities near the surface of 
formation. 
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Figure 7. - Propagation characteristics of the one-dimensional mass 
densities, ^ ■ 1.70 atm. 
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Figure 9. -One-dimensional mass flux distributions. 
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Figure 10. - Schematic formation-propagation characteristics of mass density 
functions with no vaporization. Vp >u. 
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